Manipulating Higher Partial- Wave Interatomic Interaction with an Optical Feshbach 

Resonance 
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We demonstrate a p-wave optical Feshbach resonance (OFR) using purely long-range molecular 
states of a fermionic isotope of ytterbium 1,1 Yb, following the proposition made by K. Goyal et al. 
[Phys. Rev. A 82, 062704 (2010)]. The p-wave OFR is clearly observed as a modification of a 
photoassociation rate for atomic ensembles at about 5 /jK. A scattering phase shift variation of 
5n — 0.023 rad is observed with an atom loss rate coefficient K — 7.5xl0~ 12 cm 3 /s. 



Ultracold atoms have recently been utilized as a ver- 
satile test bench for the study of many-body physics. 
In particular, the capability to tune a wide range of 
interatomic interaction with fine controllability offered 
by magnetic Feshbach resonances (MFRs) 0, 1 IS so 
powerful and has enabled numerous novel observations, 
including a Mott insulator of fermionic atoms in an 
optical lattice jij], Bose-Einstein condensation(BEC)- 
Bardeen-Cooper-Schriefer(BCS) crossover 0-0, and Efi- 
mov trimcrs[8[. For these studies, alkali atoms equipped 
with the MFR have been the main work horse. 
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An alternative approach, using an optical transition 
to artificially form a Feshbach resonance, optical Fesh- 
bach resonance (OFR) was proposed I9| and successfully 
demonstrated using alkali atoms [10h12| . These exper- 
iments were accompanied with a rather large two-body 
inelastic atom loss due to the photoassociation (PA). This 
loss can be mitigated, however, using a narrow intercom- 
bination transition in the alkaline-earth-likc atoms, as 
suggested by Ciuryfo et al. [13j . This narrow line OFR 
has been demonstrated experimentally in thermal gases 
and condensates of ytterbium(Yb) [HI, EH, which showed 
about an order of magnitude suppression of the atom loss 
as compared to the case for the alkali atoms. Later exper- 
iments using strontium atoms also showed similar advan- 
tages [l6, 17 1 . Besides the suppression of an atom loss, 
the optical manner to control the interatomic interaction, 
including the optically controlled MFR [l8|, introduces 
new possibilities in quantum gas manipulation such as 
a fast temporal manipulation and fine spatial resolution, 
which have been successfully demonstrated in Ref . fl5| . 
Arbitrary control of interatomic interactions among dif- 
ferent kinds of atom pairs in a mixture of gases would be 
also an interesting possibility. These technical advance- 
ments can broaden the spectrum of experiments that can 
be performed with quantum gases. 

In all of the previous experiments, the effect of the 
OFR has been examined only on the s-wave interaction. 
In this paper, we extend the ability of the OFR to control 
interatomic interaction of higher-partial waves 
We successfully demonstrate the p-wave OFR effect in 
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FIG. 1: (color online) Relevant molecular states and adiabatic 
potentials for the p-wave OFR experiment. Tuning and mon- 
itor lasers (A = 555.8 nm) are tuned near PA resonances at 
-355.4 MHz and -212.4 MHz from the 1 So+ 3 Pi asymptote. 
Condon radius of tuning and monitor states are R c = 5.96 
and 7.22 nm, respectively. The excited states are previously 
observed PLR molecular states arising from the hyperfine in- 
teraction [lj . A part of the wavefunction of the p-wave scat- 
tering state ^ S cat as well as the wavefunction of the state for 
the monitoring f moll are shown. Note that the energy scale 
of the ground state potential is enlarged by a factor of hun- 
dred compared with the excited state potential to show the 
centrifugal barrier of approximately 44 fiK in height. 



fermionic 171 Yb atoms in the vicinity of the purely long- 
range (PLR) molecular resonance , the use of which has 
been suggested by Goyal et al. [2fJ. An intriguing ap- 
plication of a p-wave OFR would be the study on the 
p-orbital bands of spinless fermions trapped in an optical 
lattice 2(], 21 1 . Such p-orbital physics in optical lattices 
has been discussed in Refs. [22, 23|. 



The radial part of the energy-normalized ground-state 
p-wave scattering wavefunction has an asymptotic form 
^ sca .t(r,k) = (2/z/7rft 2 fc) 1/2 sin(fcr - 7r/2 + 77) at a long 
interatomic distance r, where 77 is the scattering phase 
shift, fi is the reduced mass of the scattering atom pair, 
and k is the wave number of the atom pair, approximated 
as k = \f2fxkBT /h with T the temperature of the atomic 
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FIG. 2: (color online) (a,b) p-wave OFR spectra with only the tuning laser (solid circles) and with both the tuning and 
monitoring lasers (open circles). Number of atoms remaining after the OFR pulse (N) is plotted with respect to the tuning 
laser detuning from the 1 So+ 3 Pi asymptote (A a t om ). The data obtained with different temperatures, (a) T = 7.2 fiK and (b) 
4.5 /"iK, are shown. The dotted lines show the atom loss spectra calculated with the laser induced linewidth r tun /27r = 16.9 kHz 
and 12.3 kHz, for (a) and (b), respectively. The asymmetric atom loss spectra, indication of the OFR effect, is clearly observed 
with the monitor laser. The spectra are well reproduced in solid lines with minimal fitting parameters as discussed in the main 
text, (c) Calculated variation of r\ using the value of r tun obtained from the fit for T = 7.2 /jK. 



sample, ks the Boltzmann constant, and H the Planck 
constant divided by 2ir. The OFR effect is induced by 
the laser excitation to a molecular state, which alters the 
ground state scattering wavefunction. U sing the semi- 
analytic formalism by Bohn and Julienne 2J] , a scatter- 
ing matrix element S = cxp(2iri) for the colliding ground 
state atoms modified by the OFR laser excitation is given 
by 



S = exp(2i?7o) 



A - i (r - 7 ) /2 
A + i (r + 7) /2' 



(1) 



where 770 is the scattering phase shift associated with the 
unperturbed scattering wavefunction. 7 is the radiative 
decay rate of the excited molecular state, and A is the 
detuning of the OFR laser with respect to the molecular 
PA resonance. T is the laser-induced width given by 



Isat 



r = 77 ( 7- ) ft 7a/FC, 



(2) 



where I is the OFR laser intensity, I sa t — 0.14 mW/cm 2 
and 7 a /27r = 182 kHz are the saturation intensity and 
linewidth for the atomic 3 Pi state, respectively, fpc = 
|(V> c |d • ei|'0scat)| 2 /2(i 2 i is the Franck-Condon factor in- 
cluding the rotational correction, where d, cLa, and are 
the dipole moment operator, atomic dipole moment and 
the laser polarization, respectively. ij) e and V'scat are the 
wavefunction of the excited state and scattering state, 
including the spin and rotational degree of freedom. The 
variation of the scattering properties by the OFR effect 
can be determined from the S matrix. The loss rate co- 
efficient K and scattering phase shift 77 can be derived 



from 5-matrix as, 



(1 \ S \2) - nh r ^ 



K = 
V = arg(S)/2 



(3) 
(4) 



For the demonstration of the p-wave OFR effect, we 
tune the OFR laser near to the previously observed PLR 
molecular state arising from the hyperfine interaction 
The interatomic potential connecting to the 1 S'o+ 3 Pi 
asymptote is approximately 750 MHz in depth. Relevant 
molecular states and adiabatic potentials are denoted in 
Fig. [TJ The PLR molecular states have the inner turning 
point at ^2.5 nm and have small Franck-Condon over- 
lap with tightly bound states, possibly suppressing the 
inelastic atom loss. In order to determine the variation 
of the scattering phase shift in a systematic manner, we 
use the same technique previously demonstrated for the 
measurement of the s-wave OFR effect 0, El- Along 
with the strong OFR laser (tuning laser) which tunes 
the scattering phase shift, we apply a much weaker laser 
(monitor laser) tuned right on another PA resonance. PA 
resonances at —212.4 MHz and —355.4 MHz are associ- 
ated with T c = 3 PLR states which are optically accessi- 
ble from the p-wave scattering states, and are used for the 
monitor and tuning lasers, respectively, where T a is the 
total angular momentum of the molecular states. Con- 
don radius of the tuning and monitor states are R c = 7.22 
and 5.96 nm, respectively. 

In the previous demonstration of the s-wave OFR, the 
tuning laser varies the s-wave scattering phase shift, ef- 
fectively translating the scattering wavefunction in the 
radial direction. This translation results in the variation 
of the Franck-Condon factor for the monitor PA transi- 
tion, which enables to observe the variation of the phase 
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shift 77 in terms of the monitor PA rate variation. For the 
current study of the p-wave OFR, however, we should be 
careful due to the existence of the centrifugal barrier as 
discussed later in details. 

The experimental setup is nearly the same as our pre- 
vious works on the PA experiments and we only explain it 
briefly [25], [26|. For the preparation of cold 171 Yb atoms, 
a beam of Yb atoms is first slowed down by a Zeeman 
slower, and then loaded into a magneto-optical trap oper- 
ated with the intercombination 1 Sq- 3 Pi transition. Next, 
the atoms are transferred to a dipole force trap (532 nm) 
with a cross-beam configuration. We perform evapora- 
tive cooling by lowering the power of the dipole force trap 
lasers. A typical atom number obtained after the evapo- 
ration is 1.7 x 10 5 at 5.6 /iK. After the preparation of the 
cold atom sample with mixed spin, the tuning and moni- 
tor lasers with A = 556 nm are turned on to perform the 
OFR experiment. The two beams are overlapped and fo- 
cused at the atomic cloud with the beam waist of 70 /im. 
The tuning laser power of 525 fiW with the pulse length 
of 30 ms is used. The monitor laser frequency is fixed 
to the peak of the PA resonance at —212.4 MHz, while 
the tuning laser is scanned around the PA resonance at 
—355.4 MHz. The atom loss induced by the two lasers is 
monitored by taking the absorption image of the sample. 
We perform the measurements at a nearly zero magnetic 
field, where the Zeeman sublevels in the excited state are 
degenerate within the linewidth of the excited state. 

The observed atom loss spectra are shown in Fig. [2] 
for (a) T = 7.2 /iK and (b) 4.5 /xK. The p-wave cen- 
trifugal barrier for Yb atoms is approximately 920 kHz 
(44 /iK), and we were unable to observe clear PA spec- 
tra below 3.5 /iK. With the tuning laser only, the spectra 
are symmetric as expected for the ordinary PA resonance. 
When the monitor laser is turned on, the additional atom 
loss due to the monitor laser PA is introduced, which is 
observed for the data at off- resonance in Fig. 2(a) and 
2(b). The OFR effect manifests itself around the tuning 
laser PA resonance, where the dispersive-shaped spectra 
are observed. This behavior suggests that the Franck- 
Condon factor between the scattering wave state and the 
monitor state is altered dispersively around the tuning 
PA resonance, with the effect diminishing at a lower tem- 
perature. 

For the quantitative analysis of the OFR spectra with 
the monitor laser on, it is important to analyze the over- 
lap integral between the scattering wavefunction of the 
ground state Vfscat and that of the molecular bound state 
used for the monitoring PA. The monitor-laser- 
induced atom loss rate coefficient K mon is proportional 
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to r mon when the laser is tuned on-resonance (A n 



0) 



with the weak excitation r mon <C 7, as understood from 
Eq.[3l 

As shown in Fig. 1, the wavefunction of the molecular 
bound state $ mon is mostly inside the centrifugal barrier 
(r <7.0 nm) of the p-wave scattering state. The cal- 
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FIG. 3: (color online) Calculated wavefunction overlap inte- 
gral, ( v I'mon| , I'scat}, for the monitor state (—212.4 MHz) and 
the scattering state as a function of the phase shift 77, for 
(a) a wide range and (b) near the background phase shift 
770 = 0.127 at T = 5.6 /iK. The overlap integral is normal- 
ized with respect to the overlap integral at 770, denoted as 
(^monl^I'scat)o- Dots are the numerical results with lines show- 
ing the guide to the eye. Large overlap integral can be seen 
near 77 = 1.7 and —1.4, representing the location of the shape 
resonances. Near 770 the variation of the overlap integral is 
quite linear with respect to 77, with zero at 77' = 0.116. The 
shaded region shows the range where the variation of 77 is 
observed at a current study. 



eulated wavefunction overlap integral, (5 , m0 n| , J r scat) as 
a function of the scattering phase shift 77 is shown in 
Fig. [3]Ja). The overlap integral varies sinusoidally with 
respect to 77. The background scattering phase shift 
calculated for 171 Yb is 770 = 0.127 at T = 5.6 /iK. 
As shown in Fig. EJb), the variation of the overlap in- 
tegral is quite linear in the vicinity of 779, with zero 
overlap at 77' = 0.116. Therefore, we expect the form 
of the monitor-laser-induced atom loss rate coefficient 
Kmon = Kq{t] — r]') 2 for the current experiment. 

The data with the monitor laser is fitted using an in- 
elastic rate equation h = — 2/itotai'T 2 , where n is the atom 
density, with a PA rate coefficient if total = A' tun + K mon . 
We fit the atom loss spectra with fitting parameters r tun 
and Kq, with 7/27r = 364 kHz used for the molecular 
state radiative decay rate, and the results arc shown in 
solid lines in Fig.[Ua) and[2Jb). The calculated results fit 
the data remarkably well. Although the fit slightly over- 
estimates the atom loss near the resonance, the overall 
fitting quality including the reproducibility of the asym- 
metric spectra is satisfiable. We obtained the best fit 
with rt im /27r = 16.9 and 12.3 kHz, where the maximum 
scattering phase shift 77 = 0.150 and 0.144 radian and cor- 
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FIG. 4: Calculated variation of the scattering wavefunction 
'I' scat for various r\. A large modification of the wavefunction 
inside the centrifugal barrier (r < 7.0 nm) can be seen within 
a small variation of r\. 



responding atom loss rate coefficient if tun = 7.5 x 10~ 12 
and 5.3 x 10 -12 cm 3 /s are calculated, for T = 7.2 and 
4.5 fiK, respectively. From r tU n obtained from the fit, 
we also calculate the atom loss spectra with a loss rate 
coefficient K tota i = K tun for the data without the mon- 
itor laser and results are shown in dotted lines in the 
same figure. The loss spectra are well reproduced with 
the fitted parameter. 

The calculated variation of r\ from the obtained value 
of r tU n at T = 7.2 /iK is shown in Fig. HJc). The phase 
shift variation dr) = rj max — r]o = 0.023 radian is observed. 
On the blue side of the resonance, the OFR induced r\ 
crosses the r\ m % n = 0.116, where we expect the Franck- 
Condon factor for the monitor PA to be zero. The atom 
loss curve with the tune and monitor laser (solid line) 
tangentially connect to the atom loss with the tuning 
laser only (dotted line) at this region. From all the data 
including that for a different temperature, we calculated 
the average optical volume [20| V op t = T/2k 3 j = 29 nm 3 
at the tuning laser intensity 1 = 1 W/cm 2 for this transi- 
tion which is slightly lower than the previously calculated 
value 43.6 nm 3 at / = 1 W/cm 2 , which is the average 
value over different projections of the partial- wave angu- 
lar momentum and the excitation photon helicity |20j . 

One may wonder why the change of the monitor PA 
rate is rather strong, showing strong asymmetric atom 
loss spectra, while the phase shift variation 5rj is not sig- 
nificantly large. In Fig. 21 we show calculated scatter- 
ing wavefunctions for r\ ranging from 0.098 to 0.174 ra- 
dian. Despite the small variation in ?y, a drastic change in 
the shape of the wavefunction can be observed inside the 
centrifugal barrier (r < 7.0 nm). At rjo, the scattering 
wavefunction is away from the shape resonance condi- 
tion, and the penetration inside the centrifugal barrier is 
small and thus the amplitude of the wavefunction inside 



the centrifugal barrier is small. When the OFR causes 
a change of the phase shift rj, even if it is not so large, 
a large modification of the wavefunction can take place 
inside to alter the Franck-Condon factor with respect to 
$ mon . By choosing the molecular state which is located 
well inside the centrifugal barrier, the monitor PA signal 
amplifies the small shift of rj, resulting in a large spectral 
modification. 

In conclusion, we successfully observed the p-wave 
OFR effect as a modification of a PA rate for atomic en- 
sembles at about 5 /iK. We have utilized the PRL molec- 
ular states for the tuning and monitoring of the p-w&ve 
OFR effect. The analysis of the data shows that the 
scattering phase shift variation of Srj = 0.023 radian. 

In order to realize a larger phase shift 77 at a low tem- 
perature, it is beneficial to use an atomic species which 
has a p-w&ve shape resonance, which results in an en- 
hancement of r tU n through the large Frack-Condon fac- 
tor. The 1 Yb homonuclear pair has a p-wave shape 
resonance at the collision energy of tens of micro Kelvin 
firjl ] , and thus it would be a good candidate for an efficient 
p-wavc OFR. A future investigation will also include a 
use of ultra-narrow transitions to 3 Po,2 states in alkaline- 
earth-like atoms to further suppress the atom loss and 
heating of the sample. 
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